Full time nonexponential decay in double-barrier quantum structures 
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We examine an analytical expression for the survival probability for the time evolution of quantum 
decay to discuss a regime where quantum decay is nonexponential at all times. We find that the 
interference between the exponential and nonexponential terms of the survival amplitude modifies 
the usual exponential decay regime in systems where the ratio of the resonance energy to the decay 
width, is less than 0.3. We suggest that such regime could be observed in semiconductor double- 
barrier resonant quantum structures with appropriate parameters. 
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I. INTRODUCTION 



The exponential decay law has been very successful 
in describing the time evolution of quantum decay P, |2| ■ 
However, almost 50 years ago, Khalfin Q showed that for 
quantum systems whose energy spectra is bounded from 
below, i.e., (0,cxd), which encompasses the vast majority 
of systems found in Nature, the exponential decay law 
cannot hold in the full time interval. The present com- 
monly accepted view of the time evolution of decay in- 
volves three clearly distinguishable time regimes in terms 
of the lifetime of the system HllllllllH: (a) Nonex- 
ponential decay at short times, (b) Exponential behavior 
spanning over many lifetimes at intermediate times, (c) 
Nonexponential decay as an inverse power law of time at 
long times. The experimental confirmation of the non- 
exponential behavior has remained elusive over decades. 
After years of experimental effort, dealing mainly with 
radioactive atomic nucleus |1C| . and elementary parti- 
cles the deviation from the exponential decay law 
in the short-time limit has been finally reported some 
years ago for an artificial quantum system |12|. In the 
framework of an exact single resonance decay model |l3j , 
it is illustrated that the deviation at long times depends 
on the value of the ratio of the resonance energy to 
the decay width F^, i.e., R = £r/^r 0|- As the value 
of R diminishes, from very large values up to values of 
the order of unity, the long time deviation from expo- 
nential decay occurs earlier as a function of the lifetime, 
T = h/Tr, of the corresponding system. In a recent work, 
Jittoh et al. have shown that for values of R still 

smaller i.e., less unity, there exists a novel regime where 
the decay is nonexponential at all times. These authors 
left the discussion of full time nonexponential decay in 
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actual physical systems for future work. 

In this work we consider an analytical expression for 
the time evolution of decay for finite range potentials to 
discuss further the regime of nonexponential decay in the 
full time interval. We show that the absence of the expo- 
nential period in decay is due to the interference between 
the exponential and nonexponential contributions to de- 
cay. It is also suggested that one-dimensional semicon- 
ductor double-barrier resonant quantum structures may 
be suitable systems to verify experimentally that behav- 
ior. 

Section II presents the formalism. Section III deals 
with the calculations and its discussion, and Section IV 
gives the concluding remarks of this work. 



II. FORMALISM 

Let us therefore consider the decay of an arbitrary state 
^^{x, i = 0), initially confined at t = 0, along the internal 
region < x < i of a one- dimensional potential V{x) 
that vanishes beyond a distance, i.e. V{x) = for x < 
and X > L. The solution ^->{x, t) at time t > may 
be expressed in terms of the retarded Greens function 
g{x,x';t) of the problem as. 



tp{x,t)= / g{x,x';t)^p{x' ,0)dx' . 



(1) 



The survival amplitude A(t), that provides the probabil- 
ity amplitude that the evolved function ip(x,t) at time t 
remains in the initial state tp{^, 0) is defined as. 



Ait) 



ip*{x, O)-0(a;, t) dx, 



(2) 



and consequently the survival probability reads S(t) = 
|^(i)p. A convenient approach to solve Eq. ||2Jl is by 
Laplace Transforming g{x,x';t) into the momentum k- 
space to exploit the analytical properties of the outgoing 
Green function G~^{x,x';k) of the problem |0. Here 
we follow and generalize to one dimension the approach 
developed by Garci'a-Calderon in three dimensions Il7|. 
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The essential point of our approach is that the full out- 
going Green's function G^{x,x']k) of the problem may 
be written as an expansion involving its complex poles 
{fc„} and residues, the resonance functions 'l8'|. 
We restrict the discussion to potentials that do not hold 
bound nor antibound states. In general, the complex 
poles kn = an - i/3„ with (a„,/3„) > 0, (n = 1,2,...,), 
are simple and are distributed al ong the lower- half of the 
A;-plane in a well known manner jig. From time-reversal 
considerations, those seated on the third quadrant, fc_„, 
are related to those on the fourth, A:„, by fc_„ = — fc*. 
Analogously, the residues fulfil u_„ = u* . The complex 
energy poles En = £„ — zr„/2 may be written in terms 
of kn as En = h^kn/^m and hence £„ = ^^(a^ — /^D/^™ 
and r„ = ?i^(4a„/3„)/2m, with m the mass of the par- 
ticle. Since the energy of the decaying particle, e„, is 
necessarily positive, the poles of the system must be 
the so called proper resonance poles i.e., poles satisfy- 
ing a„ > /?„. Note that this implies that i? > 0. As 
a result of the above considerations the survival ampli- 
tude may be expressed as a sum over exponential and 
nonexponential contributions, the latter being in general 
relevant at very short and long times compared with the 
lifetime. Hence we write llTi. 



A{t) = 



{CnCn 



^-ihklt/2m 



n=l 



[CnCnM{~yn) ~ (C„C„)*M(y_„)]}, (3) 
where the function M{yq) is defined as. 



27r 



dk = ]-w{iyq), (4) 



where yq = — ex-p{—i'K /4:){h/2my/^kqt^^^ , with q = ±n, 
and the function w{z) = exp(— z^)erfc(— iz) is a well 
known function Proper resonance poles fulfil, n/2 < 
arg yn < 3n/A. The coefficients C„ and Cn in Eq. ^ are 
given by, 

Cn= I {x,0)Un{x) dx; Cn = / tp* {x , 0) Un{x) dx . 

Jo Jo 

(5) 

The above coefficients obey relationships that are similar 
to those in 3 dimensions 



eyr 

m 



Re J2 =1' E 



Cri Cn 



0. (6) 



The resonant functions Un{x), necessary to calculate the 
coefficients given by Eq. , satisfy the Schrodinger equa- 
tion of the problem with complex eigenvalues kf-^^. They 
obey outgoing boundary conditions at a; = and x ~ L, 
given respectively by, [dun{x) / dx]x=o = — jfc„tt„(0), and 
[dun{x)/dx]x=L = iknUn{L). Alternatively, the reso- 
nance functions can also be obtained from the residues 
at the complex poles of G'^ {x,x' \k) This yields a 



normalization condition that differs slightly from that in 
3 dimensions, namely. 



u^{x)dx + i 



ul{0) + ul{L) 



(7) 



The set of {fc„}'s and the corresponding {Mri}'s that fol- 
low from the solution of the above complex eigenvalue 
problem, may be obtained by well known methods|l8j. 

The long time behavior of Eq. (O, i.e., much larger 
than the lifetime r = ^/^i, rests only on the M- 
functions. At long times they behave as M{yq) « 
-a/(fc,<i/2) - b/{k^t^/^) + ...„ with q ± n, and the con- 
stants a = i/[2(7ri)i/2] and b = l/[4(7ri)i/2]. Substitution 
of the above expansion into Eq. Q gives that the factor 
multiplying t~'^/^ is proportional to the term given pre- 
cisely by the expression on the right in Eq. © , and hence 
the t^^'^ contribution vanishes exactly. This leads to the 
well known long time behavior of A{t) as t~^^'^. 

We shall be concerned here in situations where the 
initial state i'ix, 0) overlaps strongly with the lowest 
energy resonant state ui{x) of the system. In such a 
case it follows from the first expression in Eq. (|SJ), that 
Re(CiCi) is the dominant contribution. Since the decay- 
ing widths and resonance energies satisfy, respectively, 
that < Ti < r2 < and ei < £2 < it fol- 
lows that the higher resonance contributions decay much 
faster and may be neglected. This simplifies our descrip- 
tion of decay because it allows to deal with the single 
term approximation of the survival amplitude A(t) given 
by Eq. This approximation also demands to make 
sure that the correct long-time behavior of the survival 
amplitude is preserved, namely A{t) ^ t~^^'^. This re- 
quires to remove the t~^^^ contribution in the M's since 
it cancels out exactly in Eq. ^ As a consequence, 
the single term approximation of the survival amplitude 
may be written as, 

A{t)^A{t),.,p+Ait)non, (8) 

where A{t)exp is, 

A{t),xp = CiC'ie-^^'^'*/^'" (9) 

and, 

A{t)non = CiCiM{^yi) ~ {CiCiYM{y-{), (10) 

where the Ai' s denote M functions where the long time 
contribution that goes as has been subtracted to 

obtain the correct long time behavior as i^^^"^ . Hence 
the survival probability becomes 

S{t) = S(i),xp + S{t)non + S(t\nU (H) 

where S{t)expi S{t)non, and S{t)int, refer respectively, to 
the exponential, nonexponential and interference contri- 
butions to the survival probability, namely, 



(12) 
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S{t)non = \A{t) 

non I 



(13) 



and, 



Sit),nt = 2\CiCi\cos{sit/h + rj + ^{t))e-^''/^'^\A{tUnl 

(14) 

where rj — arg(CiCi) and (/)(t) = arg(A(i)„on)- long 
times, A{t)non may be written as an asymptotic expan- 
sion, that we denote by ^(i)lon' whose leading term 
reads p^ . 



-/4 ^2™^^/^J^|c^|^^ ^^^^ 



A:3 J t3/2- 



Consequently, at long times, S{t)non and S{t)int, given 
respectively by Eqs. 1)13(1 and ((14|l . may be written in 
obvious notation as. 



(16) 



and, 
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FIG. 1: 'online color'. The Survival probability S{t) (solid 
line) as a function of time for a double-barrier system with 
R as indicated and ri — 16321.9 fs. The purely exponential 
behavior of S{t) (dashed line) and the long-time asymptotic 
behavior as (dotted line), are included for comparison. 
See text. 



cos{e,t/h + 7j + 7:/A)e-^^'/^^^. (17) 



III. CALCULATIONS AND DISCUSSION 

In order to study systematically the behavior of the 
survival probability with time, in addition to the reso- 
nance parameters one needs to specify the initial state. 
It is shown below, however, that if the one-term ap- 
proximation holds, then the specific form of the initial 
state is not essential to determine the behavior with 
time of the survival probability. Hence we choose for 
ipix, 0) a very simple model, namely, the box model state 
i/'(a;, 0) = (2/w)^/^ sin7r(a: — b)/'w for b < x < b + w and 
zero, otherwise, where b and w stand respectively, for the 
barrier and well widths. The corresponding box momen- 
tum k = tt/w is closer to the real part of the resonant 
momentum ai, than to any other a's of the system. 

In Fig. n wc have used a set of potential parameters 
typical of AlAs-GaAs-AlAs double-barrier heterostruc- 
tures as in the cases considered by Sakaki and co-workers 
[20], who verified experimentally that electrons in suffi- 
ciently thin symmetric double-barrier resonant structures 
decay proceeds according to the exponential decay law. 
The potential parameters are: barrier widths b = 2.5 nm, 
well width w = 6.2 nm, and barrier heights V ~ 1.36 cV. 
In the calculations the electron effective mass is m = 
0.067me, where rrie is the bare electron mass. The res- 
onance parameters of the system are: resonance energy 
ei = 0.09959 eV, resonance width Ti = 4.0325 x 10"^ 
eV. Hence R = ei/Ti = 2469.78, much larger than unity. 
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FIG. 2: 'online color'. The Survival probability S{t) (solid 
line), as a function of time in lifetime units for a double- 
barrier system exhibiting the transition from exponential to 
nonexponential behavior. Here ri = 10.69 fs. The purely 
exponential behavior of S{t) (dashed line) and the long-time 
asymptotic behavior (dotted line), are included for compari- 
son. See text. 



and the lifetime is ri = h/Ti = 16321.9 fs. The sur- 
vival probability S{t) (solid line) is calculated using Eq. 
((TTIl . For comparison. Fig. ^ exhibits also S'(i)^o„ (dot 
line), given by Eq. 116|) . and the purely exponential con- 
tribution exp(— Fit/?!) (dashed line), i.e., Eq. (|12|l with 
CiCi = 1. We see that exponential decay law stands for 
many lifetimes. The long time nonexponential contribu- 
tion becomes relevant only after 60 lifetimes when the 
value of S(t) is extremely small, most possibly beyond 
experimental verification. 
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The oscillations of the survival probability around the 
exponential-nonexponential transition in Fig. ^ are 
caused by the cosine factor appearing in the interference 
term Sint given by Eq. (|14|l . In the long time limit, i.e., 
Eq. (|17|l . the cosine factor may be expressed in terms of 
the parameters R and r as cos(i?T+?7), and it explains the 
extremely large frequency of oscillations of S{t)) around 
the exponential-nonexponential transition for i? >> 1. 

By varying the potential parameters, one obtains also 
that the onset of the exponential-nonexponential transi- 
tion in lifetime units, depends on the value R — ei/Ti. 
In fact, it occurs earlier as R diminishes. This also occurs 
in the case of the exact single resonance decay formula, 
whose only input is the value of R 13] . One sees, from 
the above cosine factor, that the frequency of oscillations 
diminishes also as R becomes smaller. This is interest- 
ing, because it may allow to design structures with ap- 
propriate parameters to exhibit noncxponential behavior 
in ranges more adequate for experiment. 

As an example of this. Fig. [5] exhibits a plot of the 
survival probability S{t) (solid line) in a case where 
R — 0.91. The potential parameters of the double-barrier 
structure are [2l|, barrier widths & = 1.0 nm, well width 
w — 5.0 nm, and barrier heights V = 0.23 eV which give: 
£1 = 0.05639 eV, Fi = 0.06151 eV and n = 10.69 fs. The 
nonexponential behavior is set now around 15 lifetimes 
and the value of S{t) is order of magnitudes larger than in 
the preceding case. Again, for comparison, S{t)'^^^^ (dot 
line) and the purely exponential exp{—iTit/h) (dashed 
line) are plotted. Note also, since R = 0.91, of the or- 
der of unity, that the frequency of oscillations around 
the exponential-nonexponential transition is much more 
reduced than in the previous case. 

Now, it turns out that by considering systems with still 
smaller values of R, leads to the regime where that decay 
proceeds entirely in a nonexponential fashion. We have 
found that this occurs for values R ^ 0.3 Figure 
|31 illustrates an example of this regime. The potential 
parameters for the barrier widths and barrier heights re- 
main the same as in the previous case, 6 = 1.0 nm and 
V = 0.23 eV, but the well width takes now the value 
w = 1.5 nm. Note that since each monolayer of semicon- 
ductor material has a thickness of about 0.25 nm [23| . 
each barrier and the well involve, respectively, 4 and 
6 monolayers. For this system the resonance parame- 
ters are: ei = 0.07025 eV and Fi = 0.40075 eV. Also 
R — 0.1753 and ti = 1.64 fs. Here the real part of the 
complex pole, ai = 0.0491 is still larger than the corre- 
sponding imaginary part /? = .03532, which means that 
the pole is proper and provides an exponentially decaying 
contribution i.e., Eq. (|12() . Although the resonance width 
is much broader than the resonance energy, the double- 
barrier system is still able to trap the particle. One sees 
that the survival probability S{t) (solid line) exhibits a 
behavior that departs from the purely exponential behav- 
ior exjp{—iTit/h) (dashed line) along the full time span. 
Note that the long time regime, S{t)non) (dot line), be- 
comes the dominant contribution only after 15 lifetimes. 
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FIG. 3: 'online color'. Full nonexponential behavior of the 
survival probability S(t) (solid line), for a double-barrier sys- 
tem as a function of time. Here ri = 1.64 fs. The purely ex- 
ponential behavior for S{t) (dashed line), and the long time 
asymptotic behavior of S{t) (dotted line), are included for 
comparison. See text. 



Hence, one may ask what originates previously the devi- 
ation of S{t) from the exponential behavior. The answer 
follows by inspection of the expression for S{t) given by 
Eq. Hll|) : the deviation from exponential behavior is due 
to the interference term S{t)int. This is illustrated in 
Fig. 0] The interference term S{t)int (dashed-dot line), 
given by Eq. , adds up a negative contribution to the 
exponential decaying contribution S{t)exp (dashed line) 
given by Eq. (|12() . to yield a nonexponential behavior of 
the survival probability S{t) (solid line) in a time span 
that for larger values of R is usually dominated by the 
exponential term. 

As pointed out before, in the above calculations we 
have considered as initial state a box solution of well 
width w. We have examined different choices of the 
initial state to see how the nonexponential behavior of 
the survival probability is affected. We have found some 
quantitative differences but the nonexponential behavior 
remains unaffected. The differences arise from the dis- 
tinct values of the expansion coefficients as illustrated 
in Fig. O for i? = 0.1753, where we compare: (a) the 
box initial solution (solid line) [Re(CiC'i) = 0.611], (b) 
the analytical exact single level resonance formulail3i| 
(dash-dotted line) [Re(CiC'i) — 1.0], and, (c) the case 
where the initial state is the resonance function ui(x) 
along the internal region of the structure (dotted line) 
[Re(CiC'i) = 2.070]. 

We believe that a possible way to test our results for 
the full time nonexponential behavior of quantum decay, 
is by means of an experimental setup analogous to that 
used by Sakaki et al. [2^ where a laser is used to create 
electron-heavy-hole pairs in the quantum well of the dou- 
ble barrier. For thin barriers, as in the example discussed 
here, these authors showed that the decay process is dom- 
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values of R may be designed in other artificial quantum 
structures as in the decay of trapped atoms by lasers ^3 ■ 
On completing this work it came to our notice a very 
recent work by Rothe et. al. , where it is reported the 
long-awaited experimental verification of the deviation of 
the exponential decay law at long times. This has been 
achieved by measuring luminescence decays of dissolved 
organic materials. A distinctive feature of this work is 
that the small value of R is induced by a local solvent 
environment. In this respect this work differs from our 
approach which refers to the decay of an isolated sys- 
tem. It is to be expected that this experimental work 
will stimulate further research in this area. 



IV. CONCLUDING REMARKS 



FIG. 4: 'online color'. Survival probability S{t) = S{t)exp + 
S{t)„on + S{t)int (solid line) for a few lifetimes to show 
that S{t)int (dashed-dot line) causes the deviation from the 
exponential contribution S{t)e.xp (dashed line) in the pre- 
asymptotic regime for the same case of Fig. See text. 
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FIG. 5: 'online color'. The survival probability S{t) for dif- 
ferent initial states: infinite box (solid line), exact single res- 
onance (dot-dashed line) and resonant function u\ (x) (dotted 
line) . For comparison, the purely exponential behavior is also 
plotted (dashed line). See text. 



inated by tunneling escape compared with the compet- 
ing radiative recombination process. The decay rate of 
electrons is then measured indirectly by analyzing the 
time-resolved photoluminescence. What is relevant here 
is that the value of i? = ei/Ti < 0.3. Clearly these 



In summary, we have found that the full time non- 
exponential behavior of the survival probability may be 
also characterized by three regimes: (a) A first regime, 
encompassing a small fraction of the lifetime of the sys- 
tem, that is dominated by the short-time behavior and 
the high resonance contributions to the survival probabil- 
ity; (b) A second regime, dominated by the interference 
contribution between the exponential and the nonexpo- 
nential terms to the survival probability; (c) A third 
regime that is dominated by the long time asymptotic 
nonexponential contribution to decay. In fact, (a) and (c) 
are regimes that are present in general in any decaying 
system. The nonexponential behavior of decay in stage 
(b) appears in systems with a small value of the param- 
eter R in the range < i? ^ 0.3. For larger values of R 
this regime corresponds to the usual exponentially decay- 
ing behavior. Our approach possesses a general charac- 
ter for decay in quantum systems, and therefore, it may 
be applied to study the transition from exponential to 
nonexponential decay, and in particular the purely non- 
exponential regime, in other suitable designed artificial 
quantum structures. 
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